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Abstract 



o 
o 

Wc consider the explosion problem in an incompressible flow introduced in [5] . We use a novel 
jjp _ estimate for elliptic advection-diffusion problems to show that the explosion threshold 
obeys a positive lower bound which is uniform in the advecting flow. We also identify the flows for 
which the explosion threshold tends to infinity as their amplitude grows and obtain an effective 
^vq \ description of the explosion threshold in the strong flow asymptotics in a two-dimensional one-cell 

flow. 



^ ■ 1 Introduction 

^ ; 1.1 The explosion problem 

The explosion problem concerns existence and regularity of positive solutions of nonlinear elliptic 
equations of the form 

-Atf, = \g(<f>), (1.1) 

in a domain Q C M. n with the Dirichlet boundary conditions: 4> = on the boundary d^l. The 
\ nonlinearity g(4>) is convex and increasing with g(0) > and 

^ f°° ds . . 

I> ■ / — -<+oo. (1.2) 



O 1 Jo 9 ^ 

On ! 

Two typical examples to keep in mind are g(s) = e s and g(s) = (1 + s) m with m > 1. The positive 
parameter A > measures the non-dimensional strength of the nonlinearity. It has been shown in 
the pioneering works of Keener and H. Keller [24] . Joseph and Lundgren [20], and Crandall and 
Rabinowitz [12] that there exists a critical threshold A* > so that (jl.ip admits positive solutions 
for < A < A*, while no positive solutions exist for A > A*. The regularity of solutions at A = A* is 
a delicate issue: the linearized problem was studied by Brezis and Vazquez in [7] in great detail. In 
particular, when the domain is a ball, and for the exponential and power nonlinearities mentioned 
above, the solutions at the critical value A* are uniformly bounded in dimensions less or equal to 
N = 9 and N = 10, respectively, while in higher dimensions they are unbounded. For more general 
nonlinearities g(s) and domains 0, regularity of solutions at A = A* in dimensions N = 2, 3 has been 
established by Ncdcv [30J, and more recently in dimension N = 4 by Cabre [8]. 
In the present paper we consider the non-selfadjoint elliptic problem 

—A(f> + u ■ Vcf) = Xg(4>) in $7, (1.3) 
6 = on dQ, 
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with a prescribed incompressible flow u(x) so that V-u = 0. Extension of the aforementioned results 
to the case when a flow is present is a natural question in the context of the original motivation for 
the study of (jl.ip as the explosion problem [151 EH EE] . Of particular interest is to understand how 
the presence of an underlying flow and its features affect the explosion limit. 

Existence of a critical explosion threshold A* (it) can be established as a straightforward gener- 
alization of existing methods. We are mostly interested here in the qualitative dependence of A* (it) 
on the flow u — whether a flow may raise or lower the explosion threshold, and in the asymptotic 
behavior of A* (it) in the limit of a strong flow. Intuitively, a flow improves mixing and interaction 
with the boundary - hence one may expect that an incompressible flow would always raise the ex- 
plosion threshold. Somewhat surprisingly, this was shown not necessarily to be the case in [5]. More 
precisely, Berestycki, Kagan, Joulin and Sivashinsky have considered in [5], the problem (|1 .3j) for a 
two-dimensional cellular flow and observed numerically that while the explosion threshold increases 
for flows oscillating on a small scale, it may actually decrease if the flow has large scale variations. 
This is because such flows may promote creation of hot spots where the explosion would happen 
faster than without any flow. The authors of [5] have also presented a formal asymptotic analysis 
and found an effective problem in the limit of the large flow amplitude. The fully nonlinear problem 
when the flow itself satisfies a Navier-Stokes type equation coupled to the explosion problem for 
temperature has been studied in [21 [221 [29] using numerics and formal asymptotics. Recently, some 
rigorous results for the behavior of the solutions to the coupled system in the regime of a strong 
gravity have been obtained in [11]. Here we derive several qualitative properties of the explosion 
threshold \*(u) in terms of the geometry and the amplitude of the flow u. 

1.2 The main results 

In the following we always assume that O is a smooth bounded domain in M n , u(x) is a C 1 (fi) 
divergence- free flow (V • u = in f2). Our first proposition establishes the direct analog of the 
classical results for (jl.ip and allows us to define the critical parameter X*(u). 

Proposition 1.1 There exists A*(n) G (0,oo) such that (i) for every < A < X*(u) the problem 
\1.3\) has a unique positive classical solution 4>\{x) such that the principal eigenvalue k\ of the 
linearized operator Mip = —Aip + u ■ Vif) — \g'(<j)\)7p is positive; (ii) if U.3\) admits another non- 
negative solution v{x) then v{x) > 4>\; (Hi) the function <f)\{x) is increasing in A; (iv) there exists 
no classical solution of i f 1.3\) for A > A* (it) . 

The proof of this result is very close to that in [12] - we present it below both for the convenience 
of the reader and since we will use some of the intermediate steps in what follows. Another reason 
to discuss the proofs for u ^ is that some of the basic results in the self-adjoint case u = rely on 
the variational characterization of the principal eigenvalue of the linearized operator M which we 
do not have when the flow is present. 

The next theorem shows that the possible creation of hot spots cannot drop the explosion thresh- 
old arbitrarily close to zero, no matter what the incompressible flow u(x) is. 

Theorem 1.2 For any domain f2 and nonlinearity g((j)) there exists Ao > so that the critical 
threshold A* (it) for hi. Sty satisfies A* > Ao > for all incompressible flows u(x) in fL The constant 
Ao depends on f2 and the function g. 

This result does not hold without the restriction that the flow u(x) is incompressible - we describe in 
Section T2.3I examples of flows for which A* may be as small as one wishes. The proof of Theorem 11.21 
involves the following uniform LP — L°° bound for solutions of the Dirichlet problem for elliptic 
diflusion-advection problems with the constant independent of the incompressible flow. 
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Lemma 1.3 Let the flow u(x) be divergence-free and let q(x) be the solution of the elliptic problem 

-Aq + U -Vq = f(x) in $7, (1.4) 
q = on d£l, 

with f(x) G L P (Q), p > n/2. There exists a constant C(£l,n,p) > which depends on p and the 
domain f2 but not on the flow u(x), so that \\q\\L°°(n) < C||/||z,p(fi)- 

These results raise several interesting open questions. First, can one identify the optimal constant C 
in (|1.4j) ? More importantly, we would like to pose as an open problem to know whether in all domains 
with smooth boundaries a flow realizing the best constant in Lemma 11.31 exists. The same question 
pertains to the smallest possible explosion threshold in Theorem II .21 does the flow minimizing A* (it) 
over all incompressible flows exist in any domain 0? If so, what are its geometric characteristics? 
Numerical simulations in [5] indicate that a natural guess that u = turns out not to be correct for 
all domains as an incompressible flow may create additional hot spots. More precisely, it has been 
numerically computed in [5] that in a very long rectangle the explosion threshold corresponding to 
a cellular flow with a certain finite positive amplitude is smaller than that corresponding to u = 0. 
However, it is not clear whether in the situation when u = is not a minimizer of A* (it), such a 
minimizer exists at all, or if minimizing flows do not exist. When it exists, how is it determined? 

Let us now fix a flow profile u(x) and consider the explosion problem (jl.3p with a strong flow 
Au(x), with a large flow amplitude A S> 1: 

- + Au ■ Vcj) = \g(4>) in n (1.5) 
4> = on dU. 

We are interested in the behavior of the explosion threshold A* (.A) for (jl.5H in the limit A — > +oo. 
Let us recall that a function ip £ i/ 1 (J7) is a first integral of u if u ■ Vvj = a.e. in O. 

Theorem 1.4 We have X*(A) — ► +oo as A — > +oo if and only if u has no non-zero first integrals 
in f#(fi). 

This theorem provides a sharp characterization of the flows capable of preventing an explosion for 
an arbitrary A > provided that the flow amplitude A is sufficiently large. The proof uses the ideas 
from [1] and [10] together with some techniques of [Bj. Not surprisingly, the explosion threshold 
tends to infinity as A — > +oo under the same assumptions as the principal Dirichlet eigenvalue of 
the operator —A + Au ■ V (see [3]), as both quantities measure the effectiveness of the enhancement 
of the boundary cooling due to the flow. 

Finally, in Section 13.31 we consider the effective problem for (jl.5p in the limit A — > +oo for 
the class of two-dimensional cellular flows. We show that in this limit, the various cells of the flow 
"do not talk to each other" . The main result of that section is Theorem 13.41 In particular, when 
A — > +oo, the explosion threshold A* (A) is close to the explosion threshold on the "largest" cell in 
Q. Moreover, the explosion threshold for each of the individual flow cells in the limit A — > +oo has 
an asymptotic description in terms of the Freidlin problem. We recall that the fast flow asymptotics 
for the parabolic reaction-diffusion equations in flows without cells has been treated by M. Freidlin 
in [16], and our results for the elliptic explosion problem in the special case of one cell flows are 
what one would expect formally from [TBJ. The most interesting and delicate new ingredient is the 
independent behavior of the solutions in various cells when A — > +00. We mention that unlike 
in [16], our proofs are not probabilistic in nature. Actually, as a by-product of the present paper, 
one can use our arguments to recover some of the results of |16j by analytic techniques. 
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Another natural variable coefficients extension of the classical results for (11. ip is to allow the 
nonlinearity to be spatially dependent - work in this direction has been recently done in [18\ 119] . 
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2 Existence and basic properties of A* 
2.1 The critical parameter 

We begin with the proof of Proposition 11.11 which is well known for u = 0. We do not assume in 
this section that the flow u{x) is incompressible. As we have mentioned, the proof is very close to 
that in [12], with some minor modifications. Let ni[u] and rj(x) be the principal eigenvalue and the 
normalized positive eigenfunction of the adjoint problem 

— A77 — V • (urj) = fJ-i[u]r] in Q (2-1) 
7] = on d£l. 

Note that fii[u] > as the operator —A + u ■ V has no zero order term. 

Lemma 2.1 The problem U.3\) admits no non-negative classical solutions for A > m[u]/g'(0). 

Proof. Since the function g(s) is convex and g(0) > we have g(s) > g'(0)s. Therefore, any 
classical solution (f>x > of (jl.3p satisfies 

-A<(> x + u- V</> A > A(/(O)0 A in Q, (2.2) 
<j)\ = on dQ. 

Multiplying (]2.2[) by r] and integrating by parts we conclude that 

V(0) Jr](f> x < J (f)\[-Ar]-u- V77] = m[u] J # A . 

It follows that for a positive solution of (jl.3p to exist we must have fi\[u] > g'(0)X and thus no 
non-negative solution of (|1.3p exists if A > fii[u]/g' (0). □ 

Next, we show that for a sufficiently small A > a positive solution of (jl.3p exists. Let t{x) be 
the expected value of the exit-time, solution of 

-At + u- Vr = 1 in Q, (2.3) 
t = on dQ, 



and let 



6> n = maxT(x). (2.4) 



Lemma 2.2 There exists a constant C > which depends only on the nonlinearity g(s) but not on 
the flow u{x) so that problem M.Sfy admits a minimal non-negative solution <fi\ for all A < C/9 U . 

The proof is by constructing a super-solution and using it to show that a positive solution of (|1.3p 
exists. Let us recall the following fact. 
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Lemma 2.3 Assume that there exists a smooth function 4>{x) > satisfying 



-A(j> + u-V<t>> Xg(4>) inQ, (2.5) 
4> > on dQ. 

Then there exists a classical solution (f>\ of U.3\) which is minimal in the sense that for any other 
non-negative solution tp of hi. 3^ we have (f>\(x) < ip(x) for all x £ f2. 

Proof. We construct an approximating sequence 4>n{x) by setting 4>o(x) = and letting (p n +i be 
the smooth solution of 

-A(f> n+ i + u ■ V<p n+ i = Xg(4> n ) in Q, (2.6) 
</) n+ i = on dQ. 

The difference w\ := <j)\ — (f>o(= <j>i) satisfies 

-Awi + u ■ Vwi = \g(0) > in Q, (2.7) 
w\ = on dQ. 

It follows from the maximum principle that w\ > and thus > (J)q. Similarly, we have for the 
higher differences w n = 4> n — <j> n -i' 

-Aw n + u ■ Vw n = \[g((j) n -x) - g(<j} n - 2 )] > in 0, (2.8) 
w\ = on dQ. 

Then by induction we conclude that < (fr n < (t> n +i since the function g(s) is increasing. The same 
induction argument shows that (fi n (x) < 4>( x ) f° r all n > 1. Therefore, the sequence <j) n converges to 
a limit 4>\ which has to be a solution of (11. 3j) and satisfy < <j)\ < (j)(x). As the sequence <j) n does 
not depend on the choice of the super-solution (j), the limit 4>\ is a minimal solution of (|1.3|) . □ 

Proof of Lemma 12.21 Observe that for A > sufficiently small the function f(x) = 2g(0) Xt(x) 
satisfies 

-At + u ■ Vf > Xg(f) in O, (2.9) 
f = on dQ. 

Here t{x) is the solution of (|2.3[) . This is true provided that 2^(0) > g{2g{Q)Xr). As the function 
g{s) is increasing, for this inequality to hold it suffices to require that 2g(0) > g(2g(0)X9 u ). This 
condition is clearly satisfied if A < C/9 U with a constant C which depends only on the function g(s). 
Now, existence of a minimal solution to (jl.3p follows from Lemma 12.31 □ 

Recall that a solution <p\ of (|1.3p is stable if the principal eigenvalue aci(A,0a) of the linearized 
operator 

M A V> = -A^ + u ■ W - Xg'((j>\)ip 

is positive. 

Lemma 2.4 ^4ny minimal solution of il.3\) /ias Ki(A, i^a) > 0. 

Proof. Let <^>a be a minimal solution of (II. 3p and assume that the principal eigenvalue ki(A,0a) of 
the problem 

-Aip + u • VV> - A 5 '(0a)^ = «i(A, </>a)V>, (2-10) 
?/> = on 917 
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is negative. Consider the function %fi e = (fix — etfi, then we have 

-Aifi £ + u ■ Vip £ - Xg(ifi e ) = ^g{(f>\) ~ e^ff '{4>\)i> ~ £Ki((f) X )if; - Xg{(fi x - etp) 

= - £Kl (X,(fix)^-^^^>0, 

provided that e is sufficiently small and Ki(X,(fix) < 0. This contradicts minimality of (fi. Therefore, 
we have ni{X,<fix) > if (fix is a minimal solution. □ 

Lemma 2.5 Assume that <fix is a solution of hi. Sty such that Ki(X,<fix) = 0. Then no classical 
solutions of il.3\) with A > A exists. 

Proof. Assume that A > A and there exists a function (fi > such that 

-A(f> + u- X7<fi = \g(4>), 
<fi = on dQ,. 

Let us also denote by ip the positive eigenfunction of the adjoint problem 

—Aip — V • (uip) — \g'((j)\)ip = 0, (2.11) 
ifi = on (90. 

Set n = <fi\ + t(4> — (fix) with r G [0, 1]. Then convexity of g implies that 

-At? + u-Vt]- Xg(rj) = -An + u-Vn- \g((fix + r(0 - (fix)) (2.12) 
> -Ar/ + u ■ Vr/ - A(l - T)g{cfix) - Xrg(4>) = (A - X)rg(4>) > 0, 

for all r G [0,1]. Moreover, we have equality in fl2. 12f) when r = 0. Differentiating (|2,12p with 
respect to r at r = gives the following inequality for £ = (fi — (fiy. 

- AC + u • VC - Ar/(0 A )C > (A - A) 5 (<£) > 0. (2.13) 

Multiplying (|2.13p by the eigenfunction ifi of (|2.1ip and integrating we obtain 

< J ifi [-AC + u ■ VC - Xg'iMC] = J C [~N> - V • (uip) - Xg'((fix)^] = 0. 

This contradiction shows that no solution of (jl.3|) for A > A may exist if K\(X, (fix) = 0. □ 

This also finishes the proof of Proposition 11.11 The critical threshold A* (u) is the supremum of 
all A for which a stable solution of (|1.3p exists. We summarize the upper and lower bounds for A* (it) 
in Lemmas 12.11 and 12.21 as 

On ff'(0) 
We will use these bounds in the sequel. 

2.2 A uniform bound away from A* 

Uniform L°°-bounds for the functions (fix* at A = A* are difficult to obtain and will be investigated 
elsewhere [3]. However, we have the following universal estimate for A < A* which will prove useful 
later. 
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Proposition 2.6 For any 5 > there exists a constant C(5) > which depends only on 8 and 
nonlinearity g(s) but not on the domain £1 or the incompressible flow u(x) so that the minimal 
positive solution 4>\{x) of ( fi.3j) satisfies < 4>\{x) < C{5) for all A G (0, (1 — 8)X*). 

Proof. The proof is based on an idea from [6]. Fix <5 G (0, 1), let Ao < (1 — 5)X* and take any 
Ai G ((1 — (5/3) A*, A*). We denote by (fio and <p\ the corresponding classical solutions of (|1.3j) with 
A = Ao and A = Ai, respectively. 
Following [6], set 

<" ds' 



h(s) h m 

It follows from the positivity of the function g(s) and (|1.2p that h(s) is an increasing positive function 
with h(+oo) < +oo. We now define the rescaled inverse function 



h(s) 



(2.15) 



Note that, since Aq < (1 — <5)A* and Ai > (1 — <5/3)A* we have 



< ^h(s) < ^h(+oc) < -1 
Ai Ai 1 - 



(5/3 



h(+oo). 



Therefore, the function $(s) is well-defined for all s > 0, and there exists a constant K(5) which 
depends only on the parameter 5 > and the nonlinearity g(s) so that < &(s) < if(<5) for all 
s > and all A G (0, (1 - S)X*). 

In addition, as g(s) > g(0) = 1, we have 3>(s) < s and 



ti h 



Ao t / > 
Ai 



— /i (s) 
Ai 



X g(Hs)) 



Ais(s) 

Hence, as g(s) is increasing and ^(s) < s, the function $(s) is increasing, with 

< $'(s) < Ao/Ai < 1. 

Moreover, <E> is concave: 

Aos'^W^'^W-s^WVCs) A 



(2.16) 



<&"( s ) 



Ai 

A ff($(s) 
Aig 2 (s) 



9 2 (s) 
gf (*{*)) -</{*) 



Ai<? 2 (s) 



< 



because g(s) is convex, $(s) < s and < A < Ai. 

Recall that is the minimal positive classical solution to (jl.3p with A = Ai and set 4> = ^(0i)- 
Using concavity of the function <3?(s) and expression (12.160 we observe that the function <j) satisfies 
the inequality 

-A0 + u ■ V4> = -^"(0i)|V0i| 2 - Q'i^Afa + &(<j>i)(u • V0i) 
= -$"(<£i)|V<£i| 2 + X&'MiMfa) > Xx&itMgtfi) = X Q g(4>). 

Moreover, as ^(0) = the function (ft obeys the Dirichlet boundary conditions <p = on dQ,. 
Therefore, the function (f> is a super-solution for (|1.3|) with A = Ao- Employing the same iterative 
procedure as in the proof of Proposition 12.21 we may then construct a non-negative solution <f>Q of 
(jl.3p with A = Ao which is smaller than (j)(x). However, by construction we have < <f>(x) < K(5) 
and the conclusion of Proposition 12.61 holds. □ 
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2.3 A uniform lower bound for A* 

We prove Theorem 11.21 in this section. Let <j)(x) be the minimal positive solution of (|1.3p : 

-A<p + u ■ V(f> = Xg((f>) in 0, (2.17) 
4> = on dQ,. 

According to ()2.14p . in order to obtain a uniform lower bound for the explosion threshold A*, it 
suffices to bound from above 6 U , the supremum of the exit time, defined by (|2.3p and (|2.4p . That is, 
it suffices to prove that there exists a constant M > so that 

On < M, (2.18) 

for all divergence free flows u(x) in Q. The constant M should not depend on the flow u(x). This 
bound is an immediate consequence of Lemma 11.31 



Proof of Lemma 11.31 

We write q(x), the solution of (|1.4p . as 



q(x) = / ip(t,x)dt. (2.19) 



The function ijj(t,x) satisfies the parabolic initial value problem 

ipt - Aip + u ■ Vt/> = in n, (2.20) 
-0(t, x) = on 90, 
■p(0,x) = f{x) in n. 

We will now show that there exists a pair of constants C > and a > so that for any incompressible 
flow u and any solution of (I2.20p with initial data f(x) we have a uniform bound 

Cp~ at 

m,x)\<—^-\\f\\ L i { n), (2-21) 

with any r > d/2. The proof is as in [TU] with a slight modification, we present the details for the 
convenience of the reader. First, multiplying (|2.20p by ip and integrating by parts we obtain 



2dt 



I = -||VV|li- (2-22) 



Using the Poincare inequality in Q we conclude that there exists a constant a > so that 

||^ 2 )||2<e-^-* 1 )||^i)|| 2 (2.23) 

for any pair of times t2>t\> 0. On the other hand, we have, using the Poincare inequality again, 
for all Kp < 2d/{d-2): 

HWv < c||W-|| 2 . 

Next, using the Holder inequality, with l/p + l/q = I we obtain: 

l = Jw < (/ h) Vp (/ 1^ 2 - 1 ^) 17 " < c\\n l fym\\l~ l,p , 
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provided that 

2 .l) q= ( 2 .l)j^ = 2 JLll < « 



p/ \ p J p — 1 p — 1 d — 2 ' 
or, equivalently, that p > (d + 2)/4. Therefore, we have the following Nash- type inequality in fi: 

4p/(2p-l) | L/ ,|| S +2 



II 

with s = 2/(2p — 1). However, incompressibility of the flow, the Hopf lemma and the boundary 
conditions imply that ||"0(i)||l — 11/11 1- If follows that 



Going back to (|2.22p we conclude that 

d_ 

dt 

Therefore we have a bound 



2 = _ 2 ||v^||2 < _c?Jia-. (2.24) 



s+2 

s 
1 



c 

2 < -T77II/II1- 



Combining this inequality with (|2.23[) evaluated with t± = t, i 2 = 2t\ we conclude that 

n p -at 

wmh < -^Wfh, (2-25) 

with 1/s > d/4. 

Consider now the solution operator Vt '■ ipo — ► ipif)- We have shown that 

.._ „ Ce~ at 
\\Pt\\u^i? < -jpr- 

The adjoint operator to V£ is simply the solution operator corresponding to the (also incompressible) 
flow (—u). Therefore, we have the dual bound 

Ce~ at 

Wt Wl^iP < 

which in turn implies that 

.. Ce~ at 
\\P t \\ L ^L°° < 

Putting these bounds together we obtain 

CeT at n „„ 

which is (|2.2ip . The maximum principle also implies that we have a trivial bound ||-0||x,°o < - 
Interpolating between these two bounds we get the estimate 

for any e > 0. Now, for any p > n/2 we may choose e > sufficiently small so that the kernel would 
be integrable at t = 0, and (|2.19|) would imply that < C||/||lp and the constant C > is 

independent of the incompressible flow u. This finishes the proof of Lemma 11.31 and hence that of 
Theorem O □ 



VtfWoo = 11^/2^/2/ Hoc < ll^ , t/2llL2^ jLoo ||7 3 t / 2 ||ii_L2||/||l < 2/ ||/ 
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Explosion threshold in compressible flows 

As we have mentioned in the introduction, without the incompressibility constraint the explosion 
threshold may be arbitrarily small. Indeed, according to Proposition 12.11 we have an upper bound 
X*( u ) — IJ>i( u ) Therefore, to see that no uniform in the flow lower bound for A* in compressible 
flows exists, it suffices to construct flows u n (x) such that the principal eigenvalue [J>i(u n ) — ► 0, as 
n — > +00. Such example is provided by the radial flows u n {x) = 4nx, say, in two-dimensions: 

-A0 n + Anx • V0 n = n n c/) n , 4> n >0 in B(0, 1) C M 2 , (2.27) 
<j) n = on \x\ = 1. 

Then fj, n < Ce~ cn — > as n — > +oo - this can be seen either from the general theory in [171 [25] or 
by an explicit computation. Indeed, setting (j) n = e™^' ^ n we obtain a self-adjoint problem for ip n : 

-A^„ + 4n 2 \x\ 2 ^ n = (u n + 4nKVn, > in B(0, 1) C R 2 , (2.28) 
ip n = on \x\ = 1. 

Hence, fi n satisfies the variational principle 

f \V?p\ 2 + 4n 2 f \x\ 2 \iP\ 2 



fi n = —An + inf 



In addition, we have u n > 0, as follows from the maximum principle applied to (12. 27ft For a test 
function of the form ip(x) = e~ n ' x ' q(x), where < q(x) < 1, q(x) = 1 for < \x\ < 1/2 and q(x) = 
for \x\ > 3/4 we obtain by a straightforward computation 

= -« + 0(--) = S^^S? + (e -») . 4„ +0(e -»), 
with c > 0. Therefore, // n — > as n —> +oo and hence A* —> as well. 

3 The strong flow asymptotics 

In this section we consider the elliptic problem fjl .3j) when the advecting flow is strong. Accordingly, 
we introduce a large parameter A>1 and re-write (|1.3p as 

-A0 A + Au ■ V0 A = \g(<p x ) in Q, (3.1) 
4>\ = on 

We are interested in the behavior of the solution 4>(x) of (|3.ip for large A, as well as in the dependence 
of the explosion threshold A* on the amplitude A. With a slight abuse of notation we will denote 
here by A* (A) the explosion threshold of the problem (|3,ip . 

3.1 Equidistribution on the flow streamlines 

Our first result shows that, when the flow is strong, solution becomes nearly constant on the flow 
streamlines, at least in an average sense and for A away from X*(A). This is a common phenomenon 
in diffusion-advection problems: a strong flow induces stratification. 
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Proposition 3.1 Assume in that u ■ n = on the boundary d£l. Then the solution <f)\ of $1.3\) is 
nearly constant on the streamlines of u for sufficiently large A in the sense that for any 5 > there 
exists C{5) so that 

\u- V0 A | 2 < 
/or a// A < (1 - 5)X*(A). 

Proof. First, we multiply (|3.ip by u • V0a and integrate over Q. The emerging integrals over the 
boundary vanish since (f>\ = u ■ n = on d£l, and we obtain the following estimate: 

\u ■ V0 A | 2 = jj(v V0a)[A^ a + A 5 (0 A )] = jf( u - Wa)A0a = -4 j 



A ,/q <9x fc 



1 /" 50a 3itj c^a 1 f d(p\ d 2 4>\ If dcj)\ duj d(j)\ C 



2 



A Jq dxk dxk dxj A Jq dxk J dxjdxk A Jq dx^ dxk dxj A J ^ ^ 

This means that the variation along the streamlines is smaller than across. Now, we have to bound 
the L 2 -norm of V^>>. However, multiplying (|3.ip by <f>\ and integrating by parts we see that 

|V0 A | 2 = xf gi^x. 

Moreover, Proposition 12,61 shows that <f>\ satisfies a uniform bound < (f)\ < C(5) as long as 
A G (0, (1 — 5) A*). Therefore, for such A we know that 



for all A > 0. It follows that 

CX 



J | V0| 2 < CX, (3.2) 
u- V(f>\ 2 < 



A 

In that sense solution becomes uniform over the streamlines. □ 

As a consequence, for each fixed A < limsup^,^^ A* (A) we know that 



/ 



u ■ V(/>a| 2 — ► 0, as A — > +oo. 



We will improve this statement for two-dimensional cellular flows in Section! 

An interesting by-product of the estimate f|3.2[) is that there are no boundary or internal layers 
in this problem unlike in the problems with boundary forcing in cellular flows [9l [T3"| [22l [271 13"TI [33| 
33 E7J- The reason is that 4>\ is set to be constant on the boundary and the normal component of 
the flow vanishes on the boundary - hence, there is no "conflict" between the boundary data and 
uniformization along the streamlines of the flow. 



3.2 The critical parameter in the limit of a strong flow 

We prove here Theorem 11.41 We recall that the assumption that there is no Hq(Q) first integral 
H (x) such that u ■ VH = almost everywhere is equivalent to the fact that the principal Dirichlet 
eigenvalue /J-i(A) of the operator — A + A ■ V on Q tends to infinity as A — > +oo [3]. 

First, assume that Hi(A) is bounded as A — > +oo. Then the upper bound in (|2.14p for A* (A) 
implies that limsup^,^ A* (A) is finite as well. 

Next, we show that if u has no first integral in H^iQ) then A* (A) — > +oo as A — ► +oo. The 
proof is based on the lower bound for X*(A) in (|2.14|) . The next lemma is contained in [25] - we 
provide a proof in the spirit of [10] . 
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Lemma 3.2 We have 9 a — ► as A — ► +oo i/u /tos no /irsi integrals in Hq(Q). 

Proof. As in the proof of Theorem 11.21 we represent the function ta(x) using the Duhamel formula 

as 

t a (x) = / i/j(t,x)dt, 



Jo 

with the function ip(t,x) that solves the parabolic problem 

tpt - A^> + Au ■ Vip = in O, (3.3) 
V>(0,x) = 1, 
tp = on cKl 

It follows from the Theorem 5.3 of [10] that if u has no first integrals in Hq(£1) then for any to > 
we can find a flow amplitude Aq^q) > so that ||V ; ( n ^o)||L° o (n) ^ 2~ n |f2| for all A > Aq. Therefore, 
we have an upper bound 



oo /-,\n 

t a (x)\ < t \n\ ^2 U = 2t °^ for A - A °^- 

n=0 



We conclude that UtaIIl 00 ^) ^ as j4 ^ +oo. □ 
This also finishes the proof of Theorem II. 41 □ 

3.3 The explosion problem in a two-dimensional cellular flow 

We now consider the explosion problem 

-A4> A + Au ■ V(/)a = XfiM in n, (3.4) 
4>A = on <9f2, 

in a two-dimensional simply connected domain £1. The flow has the form u = (^ y ,—^f x ) with a 
stream-function ty(x,y) which we assume to be sufficiently smooth. We assume that the boundary 
of the domain O is a level set {^f = 0} which may contain finitely many saddle critical points of 
the function ^ - thus, the boundary is a union of streamlines of the flow away from the critical 
points. We also assume that inside O the flow has a cellular structure: the saddles of \f are all 
non-degenerate and are connected by the flow separatrices which divide f2 into a finite number of 
invariant regions, called the flow cells, that we will denote by Cj. The stream-function ^f(x,y) has 
only one critical point (xq,i/q) inside each of Cj, which is a non-degenerate maximum or minimum. A 
prototype example of such flow has the stream-function ^(x, y) = sin7rxsin7ry - its cells are squares 
[n, n + 1] X [m, m + 1] with integer m and n, and the domain 17 is a finite union of such squares. A 
more general flow of such type is depicted in Figure [3TT1 



The Freidlin problem 

The strong flow asymptotics for parabolic reaction-diffusion equations for two-dimensional flows with 
Morse class stream-functions has been considered in [16]. This class does not include the cellular 
flows under our consideration as we allow the stream-function to have several saddles on the level 
set = 0}. Nevertheless, the limit problem of [16] is crucial in the explosion problem in a cellular 
flow. The limit problem in [16] was formulated as a system of reaction-diffusion equations on the 
Reeb graph of the function ^f. We recall and re-derive these results below in the context of the 
explosion problem in the one-cell setting, as that is what we will need below. The single cell is also 
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Figure 3.1: A schematic description of a cellular flow. 



the situation addressed numerically in [5]. For a one-cell flow the Reeb graph is simply an interval 
[0,i?o]> where Hq is the value of $f at the critical point inside the cell C that we assume to be a 
maximum, and = 0} is the boundary of the cell. We are interested in the behavior of solutions 
and of the explosion threshold in the limit of a large flow amplitude. 

The effective Freidlin problem on the interval < h < Hq is to find a function (p(h) satisfying 



1 d 



p(h)- 



T(h) dh \ jrv " / dh / 
<£(0) = 0, 4>'{h) is bounded for < h < Hq, 



(3.5) 



with the coefficients 



T(h) = I 

Jty(x,y)=h 



dl 



<(x,y)=h J^(x,y)=h 

Under our assumptions on the stream-function, the average turnover time T(h) is bounded from 
above and below away from zero: 



p(h) 



IV^IdZ. 



(3.6) 



< T < T(h) < Ti\\nh\. 



(3.7) 



The uniform bound from below by To in (|3.7p comes from the fact that the maximum of *f?(x, y) is 
a non-degenerate critical point. The term 0(|ln/i|) for small h > appears in (|3.7p because the 
boundary may contain non-degenerate saddle points of ^ so that the turnover time blows up as 
h I 0. The coefficient p(h) is positive for h > and behaves as p(h) ~ C(Hq — h) close to h = Hq. 
In particular we have p(Hq) = (diffusivity vanishes at this point), while the drift satisfies 



p'(h) 



^(x,y)=h 



IWI 



dl < — ao, with qq > 0, 
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for h near Hq, and points away from h = Hq. Therefore, the end-point h = Hq is inaccessible for 
the diffusion process corresponding to the left side of (|3.5p . and one does not need to prescribe the 
boundary condition at h = Hq in order for (j3.5l) to be well-posed. The following proposition defines 
the explosion threshold for the effective problem. 

Proposition 3.3 There exists A* > so that a positive solution of the effective problem \3. 5\) exists 
for all < A < A* and there is no positive solution of Ii3.5\) for A > A*. 



Proof. The proof follows the same steps as in Section \2. II - the only required modification is due to 
the degeneracies at h = and h = Hq. This can be addressed using the general theory in [H] and 
|28| but in the present case the boundary value problem with a prescribed right side 

1 " <v(h) d 4)=Hh) (3.8) 



T(h) dh\ i y ' dhj 
^(0) = 0, i//(h) is bounded for < h < Hq, 

has an explicit unique solution 
where 

Jo p( s ) 

In particular, we have |-P(£)| < C\ ln(//o — 01 > so that 

mh)\ < / |/(e)|r(OP(min(/i,0)de<C||/||oo / |lne||ln(iJ- -6Me<C||/|U, 
J o J o 

and we also have 

wm - m r mmm - ^ r 1 ,nf ,df - cuiu 

so that H^Hjyi.oo < C||/||x,oo. Therefore, the mapping f(h) — > ip(h) is a compact map on C[0, Hq] 
and the Krein-Rutman theory applies to the operator in the left side of (|3.8p . We may then repeat 
the proof of Proposition 11.11 essentially verbatim and conclude that the critical threshold A for (|3.5p 
exists. □ 



The explosion threshold for strong cellular flows: the main result 

The main result of this section is the following theorem. We assume that the flow has a cellular 
structure and satisfies the assumptions outlined at the beginning of this section. Then for each cell 
Cj one may formulate the corresponding one-cell Freidlin problem (13. 51) for a function <f>j, posed now 
on an interval [0, if,], where the outer boundary of Cj is the level set = 0} and Hj is the value 
of the function ^ at the (unique) extremal point inside Cj. For the Freidlin problem the Dirichlet 
boundary condition 0(0) = is prescribed at h = 0, and the derivative <j)'j(Hj) is imposed to be 
bounded. This defines the explosion threshold A* for each cell Cj. The following theorem shows that 
in the limit of a large flow the explosion threshold for the whole domain approaches the Freidlin 
explosion threshold for the "largest" cell Cj. 
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Figure 3.2: The value of A* for various values of the Peclet number Pe = A for the cellular flow with 
the stream-function tp(x, y) = sin(2-7r((2x/3 + l) 3 /8 + 1)) sin(2-7r((y + l) 2 /4 + 1)) on the the domain 
[0, 2ir] x [0, 2ir] with four cells. Dashed-dotted line: A* for the whole domain [0, 2tt] x [0, 2tt], solid 
line - the minimum of A* of the four individual cells. 

Theorem 3.4 Let A* (^4) be the explosion threshold for ^3. 4\ ) and A* be the threshold for the afore- 
mentioned effective problem \3. 5\) posed in the cell Cj of the domain £1. Then we have 

lim A*(^4) = minA*. 
A-+oo j 

A numerical illustration of the main result of Theorem 13.41 is depicted in Figure 13.21 

The proof of Theorem 13.41 proceeds in several steps. First, we prove a stratification lemma 
for solutions of forced advection-diffusion problems in cellular flows. It shows that solution of 
the Dirichlet problem is small not only on the outer boundary but also on the whole skeleton 
of separatrices and cells " do not talk to each other" . The second step is to establish the result of 
Theorem 13.41 for domains consisting of one cell where one just has to show that for one cell the 
explosion threshold converges in the strong flow limit to that of the Freidlin problem. The last step 
is to generalize this result to a domain consisting of finitely many cells. 

3.4 Cellular flows: a stratification lemma 

We begin the proof of Theorem 13.41 with the following lemma, of an independent interest. Let T>q be 
the union of all cell boundaries (separatrices of the flow) inside f2 including the outside boundary 
d£l. First, we show that solutions of a linear problem with the homogeneous Dirichlet data on the 
outer boundary are small on T>q. 

Lemma 3.5 Let i^a{x) be the exit time from £1, solution of 

—Aip a + Au ■ VipA = 1 inn, (3.9) 
ipA = on dVt. 
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For any 5 > there exists Aq > so that for all A > Aq we have < i/ja(x) < $ f or M % G T>q. 

Intuitively, this lemma says that once a diffusive particle obeying an SDE 

olX t = Au{X t )alt + V2dW t , 

comes close to the skeleton of separatrices, somewhere inside f2, it exits the domain f2 after a short 
time. This is the phenomenon behind the effective diffusivity [HI [131 E2J EH EH E21 ESI EZ] , and front 
speed and principle eigenvalue enhancement [H [26l EH Ell E9j in cellular flows. 

Proof. The proof is in two steps. First, we show that for any 5 > there exists a small h > 
and a large Aq > so that for any A > Aq we can find hj{A) with /i < \hj\ < 2h, such that 
< iPa(x) < <V2 f° r all x on the streamline {x G Cj : 'I'(x) = inside the cell Cj. Let 

be the interior of those streamlines. In the second step we consider the water-pipe domain 

Pa = Q \ ( Uj Dj (A) j , that is, a narrow tube around the skeleton of separatrices. Using the fact 
that hj is small, we apply the maximum principle for narrow domains to conclude that the function 
ipA is smaller than 5 in all of Pa and not only on its boundary. As a consequence, Va is small also 
on the skeleton Dq. 

Step 1. We have a uniform L 2 -bound for the gradient: 

\Vip A (x)\ 2 dx < C, (3.10) 

n 

with the constant C independent of A > 0, which follows from multiplying (|3.9p by iPa{x) and 
integrating by parts, together with the uniform L°°-bound for iPa(x). < C, which follows 

from Lemma 11.31 

Now, take s G (0, /i/4) and let Fj(h,s) be the domain between the two streamlines {^(x) = 
5h/4 — s} and {^(x) = 7h/4 + s} inside the cell Cj. We multiply (|3.9|) by (u ■ Vi^a) and integrate 
over Fj-. 



dipA 9 2 tpA A 

0X m OXu 



•* F i Jf 3 m,k 

= U (u-VM(n-^ A )dl-U ^?tA^A dx <C [ | V ^(x)| 2 dx + ^. 

A JdFj A J Fj dX k OXm OX k A JdF j A 

We used (|3.1U|) in the last step. Averaging this estimate in s G (0, h/A) we conclude that 

C(h) 



Pi 



\u ■ VipA.\ dx < 



A 



where Fj = Fj(h, 0) is the domain between the streamlines {^(x) = 5h/4} and {^(x) = 7/i/4} inside 
the cell Cj. The constant C(h) may blow-up as h [ but that is not important at the moment. 

It follows that there exists a value hj(A) G (5/i/4, 7/i/4) so that along the streamline Lj(A) = 
{x G Cj : ^(x) = hj(A)} we have 

I \u-V^ A \ 2 dl<^-, 

JLj{A) A 
with a new constant C{h). Therefore, the oscillation of ipA along Lj(A) is small: 

C(h) 



< 



A 
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Hence, iPa(x) is close to a constant Mj(A) on the streamline Lj(A) when A is sufficiently large. As 
a consequence of the gradient bound (|3.10p . we have \Mj(A) — M m (A)\ < C\fh if the cells Cj and 
C m have a common piece of the boundary. As the outer boundary <9f2, where ipA{x) = 0, is also part 
of some cell boundaries, it follows that < Mj(A) < C\fh for all cells Cj. Therefore, we have 

< i/> A (x) < CVh+ ^= < i, for x G Lj(A) 

if h G (0, ho) is sufficiently small and A > Aq is large enough. 

Step 2. Now, we look at the water-pipe P A and show that solution is below 5 everywhere in 
Pa- The function ipA inside Pa satisfies < iPa(x) < S/2 + r A , where is the exit time from the 
slightly larger domain Q2/1 = ^ \ (\Jj Dj(2h)J : 

-Am + Au ■ Vr A = 1 in Q 2 h, (3.11) 
r A = on (9Q 2 /i- 

Now, as in the proof of Lemma 11.31 we conclude that there exists a constant C{h) such that 
||rA||^oo(Q 2h ) < C{h) for all A > 0. The same proof shows that C{h) — > as h — > - this 
happens because the principal eigenvalue of the Dirichlet Laplacian in tends to infinity as h — > 
while the constants K p (h) in the Poincare inequality ||V ; llLp(Q 2 h) — ^p(^)ll^^lli 2 (Q2h)' 1 < P < °°i 
satisfy iT p — ► as /i — > 0. It follows that if we take h > sufficiently small (independent of A) then 
< ta{x) < 5/2 for all A > 0. Therefore, we have < iPa(x) < <5 for all A > in and, in 
particular, < iPa{x) < 5 on T>q. The proof of Lemma 13.51 is now complete. □ 

3.5 Explosion problem in a one-cell domain 

We now consider the explosion problem 

+ Au ■ V(j) A = \f(<f> A ) in ft, (3-12) 
(f>A = on <9f2, 

in a domain f2 which consists of just one flow cell. Without loss of generality we assume that the 
single critical point (xo,yo) £ of $ inside O is a maximum and set Hq = ^(xo,yo)- Let us now 
formulate the version of Theorem 13.41 for a one-cell domain. 

Proposition 3.6 Let Q be a one-cell domain and let A* (.A) be the explosion threshold for Ii3.12\) and 
A* be the threshold for the aforementioned effective problem \3. 5\) posed on [0, Ho]. Then we have 

lim \*(A) = A*. 

A-^oo 

The proof of Proposition 13.61 is in two steps. First, passing from the problem on the cell to the 
Freidlin problem we show that the Fredlin threshold A* is not smaller than than limsup j 4^ +00 \*(A). 
Next, we establish the opposite inequality by starting with a solution to the Freidlin problem and 
constructing a super-solution for (I3.12p . The second step is quite straightforward in the case when 
the boundary of SI contains no saddles of the flow u but is somewhat more technical if <9fi contains 
such fixed points. 
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Passage from the cell to the Freidlin problem 

We first prove that 

lim sup A* (yl) < A*. (3.13) 

A— >oo 

Assume that 

A < limsup A*(yl). (3.14) 

A^oo 

We will show that then A < A* by constructing a solution to the Freidlin problem (|3,5|) as the limit of 
a sequence of problems on Q. It follows from (I3.14p that there exists 5 > and a sequence A n — > +oo 
such that A < (1 — 5)X*(A n ). Therefore, as a consequence of Proposition 12.61 the minimal positive 
solutions of 

-A<p n + A n u ■ V0„ = A/(0 n ) in fi, (3.15) 
(f> n = on 50 

are uniformly bounded in L°°(Q) n i?o(fi): 

<(/>„< C, y |V(/>„| 2 (ix < C, (3.16) 

with the constant C > independent of n. Hence, the sequence </> n converges weakly in H 1 ^) (after 
extracting a subsequence) and strongly in L 2 (Q) to a function (f>. As the functions n are uniformly 
bounded and g(cfr) is smooth, the sequence g{4> n ) converges to g{4>). 

We claim that <p depends only on the variable h = ^(x,y) and satisfies the effective Freidlin 
problem (|3.5p . The first claim follows after we divide (|3.15p by A n and let n — > +oo. This leads to 

u-V<^ = (3.17) 

in the sense of distributions. It is convenient now to introduce the curvilinear coordinates (h, 9). The 
coordinates are chosen so that h(x, y) = ^(x, y), that is, the streamlines of the flow are {h = const}, 
and the level lines of the coordinate 9 = @(x, y) are orthogonal to the flow lines: V0 • = 0. We 
normalize 9 so that < 9 < 2tt and the boundary d£l is a level set: d£l = {h = 0}. Then (|3.17p 
implies that (j) depends only on the variable h. The L°°-bound in (|3.16p implies that < (j)(h) < C. 
In addition, we have 

/[•Ho _ / i-2-k |V7\T/|2 \ 
\4> h \ 2 \vh\ 2 dx = l^l 2 (/ l —^d9jdh. 

Here J = ^ y @ x — ^x&y is the Jacobian of the coordinate change. Note that V0 = pV _Ll I / with 
some scalar function p > 0, so that 

J = p|W| 2 , |V9| = p\V^\ and dl = d9/\Ve\. 

Therefore, we have 

27r IWI 2 f 

d9 = <t \V%dl =p(h), (3.18) 



J J*S>(x,y)=h 



and thus we have a weighted if 1 -bound 

-H 



p{h)\4>h\ 2 dh < +oo, 



o 
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which follows from (|3. 16j) . and hence (f>(h) is continuous for h < Ho, as p(h) ~ C(Hq — h) for h close 
to Hq. 

Next, we re- write A3. 151) in the curvilinear coordinates: 

\v*\ 2 d 2 4>n |ve| 2 d 2 ^ n (A*) 9^, (AQ)d ( f> n dK 1 

+ A n^T = -l X 9Wn), (3.19) 



J d/i 2 J d6> 2 J dh J 80 80 J 

4> n (H , 0) = 0, (f) n (h, 0) is bounded for < h < Re- 
integrating this equation in and passing to the limit n — ► +oo we obtain the limit problem for the 
function 0: 

-a(h)4>"(h) - b{h)4>'{h) = \c(h)g(0), (3.20) 
0(ff o ) = 0, 4>{h) is bounded for < h < H , 

with 

f 2n IV*I 2 /" 27r A* f' 2n 60 

<h) = / l —^d0, b(h) = / — dO, c(h) = / -. 
Jo J Jo J Jo J 

It remains only to observe that (I3.20p is nothing but the effective problem (13. 5p . Indeed, as in t|3. 181) 
we compute that 

r 2w do r di m „. 

c(h) = -=f — = T(h), 

Jo J Jv(x,y)=h l VW l 

and 

f 2n Avp r A* 

b{h)= —d0=j> dl=p'(h). 

Jo J J^(x,y)=h l VV l 

The last equality above follows from the fact that 

p(h) = i> \V^f\dl = / A^dxdy. 

J^(x,y)=h JG h 

Here Gh = {h < ^(x,?/) < -ffo} is the interior of the streamline {^>(x,y) = h}. It follows that 
(j3.20j) is, indeed, the effective problem (13. 5| , so that </>(/?.) is a positive solution of (13. 5| . Therefore, 
in particular, we have A < A* and f|3. 13j) holds. 



Subsolution: the case with no saddles on <9f2 
We now prove that 

liminf A* (A) > A*. (3.21) 

A^oo 

Together with (|3. 13j) this will complete the proof of Proposition 13.61 This is done as follows: we 
take any Ao < A* and show that Ao < X*(A) for a sufficiently large A by constructing a bounded 
positive super-solution to (|3.12p with A = Ao- However, the singular points of ^(x, y) cause technical 
difficulties in the construction of the sub-solution. Hence, we first consider the special case when 
^{x,y) has no saddles on dQ. 

Let Ao < A* and let <f>{h) be the corresponding positive solution of (|3.5|) with some A £ (Ao, A*). 
An important observation is that there exists C < +oo so that for all (x,y) ^ (xo,yo) (that is, not 
the maximum of the stream-function ^(x,y)), we have 

\A Xty 4>(h(x, y))\ < C, (x, y) + (x , yo). (3.22) 
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To see that we write 

A x , y 4> = \V^\ 2 4>"{h(x, y)) + (A*)0'(/i(x, y)) 

and note that for h close to Hq we have |V , J / | 2 ~ (Hq — h) and 4>"{h) ~ 1/{H — ho) so that the first 
term above is uniformly bounded in (x,y) G O. 

We look for a super-solution to (|3.12p in the form 

4> = 4> + rj A , 

where t]a is smooth and bounded. Then the uniform bound (|3.22[) and similar bounds on other 
second derivatives of 4>(x, y) show that if we have 

- A^ + Au- V0> A o s(0), (x,y) + {x ,y ), (3.23) 

in and <p > on then <fi(x, y) is a weak super-solution to (|3.12p with A = Ao in the sense of [6] 
and thus Ao < X*(A). We choose tja = i]A(x,y) as the solution of 

-A77 + Au- Vr] = A X)y (f> + \g((f)), (3.24) 
r](x) = on dQ. 

Assume that we can show that 

- > 0, as A — > 00, (3.25) 

then ([333]) holds: 

-A0 + An • V<£ = -A0 - Ar/ A + An • Vr, A = \g{4>) > X g (4> + rj A ) = A g(^>), 

for A sufficiently large. Above we used the fact that Ao < A and a uniform bound for cp. Hence, 
^0 < X*(A) for A sufficiently large. This will prove Proposition 13.61 in the special case when the 
domain consists of one cell and the boundary d£l contains no saddles of the stream- function *S>. 

It remains to establish fl3.25|) . To this end consider a cut-off function x( s ) so that < x( s ) — 1 
and x( s ) = 1 f° r l s l < 1/2 and %(s) = for |s| > 1 and split 

A x>y (j) + \g(4>) = qi + q 2 , 

with 

qi = (A X J + Xg®) x ( *°-f>v)) } , q2 = {Ax J + XgW) [l - X 

The small parameter <5 > is to be chosen below. We define, accordingly, the functions m, j = 1,2 
as solutions of 

-Arjj + Au ■ Vrjj = q jf (3.26) 
rjj(x) = on 90, 

so that = i]i + ?/2- As q\ is uniformly bounded in L°°(fi), the function 771 can be bounded using 
Lemma 11.31 as 

IMloo < C8 a 

with some a £ (0, 1). We may further split the function 772 = 773 + 774, with the function 773 that 
solves the ODE along the closed streamlines: 

Au-Vrj 3 = q 2 , 773(0 = 0,h) = 0. 
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This equation is solvable because 

£ [A X)V h + Xg(h)] dl = 0, 

as this is how the Freidlin problem is obtained. The function 7/3 satisfies the estimate 

IMMn) < (3-27) 
with some function Fi(5) (which may tend to infinity as 8 J, 0). Finally, 174 satisfies 

-Ar/ 4 + Au- Vt? 4 = A7/3, (3.28) 
774(2;) = — rjs(x) on dfl. 

Once again, Lemma 11.31 together with the C 2 estimate (|3.27|) on r/3 above implies that 

II II <? CFl ^ 

Altogether we see that for any e > we can find find 5 > 0, and then find Aq so that for any A > Aq 

IMU < e. (3.29) 

This proves Proposition 13.61 in the special case when the domain Q, consists of one cell and the 
boundary dfl contains no saddles of the stream- function ^. 

Approximation on a smaller domain 

Now, we establish the claim of Proposition 13.61 for domains f2 which consist of one cell but may have 
saddles of the stream- function on the boundary dQ,. In order to avoid dealing with the singular 
points on the boundary of f2 in the construction a sub-solution we need to consider a slightly smaller 
domain Vl e = {e < ^(x,y) < Hq} C O, with e > small. The domain £l e has no saddles on <9fi and 
thus the conclusion of Proposition 13.61 holds for f2 e by what we have shown above. 
Define A* (A) as the explosion threshold for the problem in £l e : 

-A(f) E A + Au ■ V<#4 = \f{<jf A ) in ft £ (3.30) 
(jf A = on dQ e . 

We also let A* be the explosion threshold for the corresponding Freidlin problem: 

1 d 



TWdfcv^r" 9 *' (3 - 31> 

<^ e (e) = 0, 4>' £ (h) is bounded for e < h < Hq, 

with T{h) and p{h) still given by (|3.6p . As we have mentioned, since O e has no saddles of ^ on d$l £ , 
we have 

lim A* (A) = A £ , (3.32) 

for all e > 0, according to the what we have already shown above. 

It is clear from the definition of A* (A) and A* that A* (A) > A* (A) and A* > A*. The next two 
lemmas show that the passage to the limit e — > is harmless. The first statement concerns the 
Freidlin thresholds. 
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Lemma 3.7 We have lim A* = A* ; where A* and A* are the explosion thresholds of $3. 5|) and $3.31\) . 

£^0 

Proof. The proof of this lemma is rather straightforward. It is clear that A* > A* for all e > 0. On 
the other hand, given 7 G (0, 1), for A < (1 — 7)A* we can find 5 > which depends only on 7 so 
that solution of the following problem exists: 

f Pih)^ = Xg{4>) (3.33) 



T(h) dh V dh 
4>(e) = 5, (fr'(h) is bounded for e < h < Ho, 

for all e > 0. Then it is easy to verify that, for a sufficiently small e > (and 5 > fixed), solutions 
of the iteration process 



p(h)^)=Xg^ n -i) (3.34) 



1 d 
~T{h)dh \f v ~' dh 
n (O) = 0, 4>' n (h) is bounded for < h < H , 

with (fro = are increasing in n, and satisfy 4> n (h) < (5 for < /i < e and (fr n {x) < for e < h < Hq. 
Thus, the sequence ^ n (/i) converges as n — > +00 to a bounded solution of (|3.5p so that A < A. □ 
The next lemma shows that A* (A) is close to A* (A) for e small. 

Lemma 3.8 For any 7 > i/iere exists £0 > and Ao so that for all e G (0, £0) a nd A > Aq we 
have (l- 7 /4)A*(A) < A* (A). 

This lemma is sufficient to show that 

A* < liminf X*(A) (3.35) 

and thus finish the proof of Proposition 13.61 for all one-cell domains, as we have already established 
(|3.13p . In order to see that (|3.35p holds, take 7 G (0, 1) and find £0 and Aq as in Lemma [3T8l Consider 
any A < (1 — 7)A*. Then Lemma [3771 implies that there exists £1 < £0 so that A < (1 — 7/2)A* 1 . 
Now, (|3.32p implies that we can find A 1 so that A < (1 - 7/4) A^ (A) for all A > Ay. As e± < e we 
may use Lemma 13.81 to conclude that A < A* (A) for all A > Aq + A\. Therefore, (|3.35|) holds. This 
finishes the proof of Proposition 13.61 □ 



The proof of Lemma 13.81 

The proof of Lemma 13.81 is based on the iteration argument and stratification Lemma 13.51 We start 
with A < (1 — 7) A* (A) for all A > Aq, £ < £0 and construct a solution of the explosion problem on 
the whole domain by the iteration procedure. Set (fro = and let (fr n be the solution of 

-A(fr n + Au ■ V(fr n = \g((fr n -i) in O, (3.36) 
(fr n = on <9f2. 

We claim that the sequence (fr n (x) is increasing in n, pointwise in x, for each A and there exists Ao 
so that for A > Ao it has a uniformly bounded limit (fr(x) G L°°(Q) which satisfies ()3.12p . 

Pointwise monotonicity of (fr n {x) in n is standard: the difference r]i(x) = (fri(x) — (fro(x) = (fri(x) 
satisfies 

-At/i + Au ■ V171 = Ao(0) > in ft, (3.37) 
rji = on (90. 
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Hence, we have r)i(x) > in Q and (f>i(x) > 4>o{x). Let us assume that i] n (x) = <j) n (x) — <ftn-i(x) > 
in ft. As the nonlinearity g(s) is increasing in s, the function rj n+ i(x) is the solution of 

-A Vn+1 + Au ■ VVn+l = Ab(</>n) - g{<t>n-\)\ > in ft, (3.38) 
rji = on dQ,. 

Therefore, we have r/ n +i(x) > inside ft and thus the sequence (j) n (x) is increasing in n. 

We need to show that the sequence </> n (x) is uniformly bounded from above. We recall that 
A < (1 — 7) A* (A) and thus the minimal positive solution of (|3.30p satisfies < <fA < C(7) with the 
constant C(~f) independent e and A. Hence, we may find 5 > which depends only on 7 > but 
not on e or A so that solution of the following problem exists: 

-A( + Au-V( = \g(() in ft e , (3.39) 
C = 5 on dfl £ . 

To see that such 5 > exists, let ip(x) satisfy (|3.30p with A' = (1 — 7/2)A*(A) and set r(x) = ip(x) + 5, 
then r(x) satisfies 

- Ar + Au ■ Vr = \'g(r - 6) in ft £ (3.40) 
r = 5 on <9ft e . 

It is a super-solution for (|3.39[) if we ensure that \'g{r — 5) > Xg(r), or, equivalently, 5 G (0, 1) is 
taken so small that we have for all x 6 ft: 

>^^. (3.41) 



<7(#c)+<5) " 1-7/2 ■ 

The function ^(x) obeys a uniform bound ||'!/'I|l oo (c) — -^(7) with if (7) independent of e and A 
Let M = sup < s <^-( 7 ) + i g'(s), then (|3.41j) is guaranteed if we have 

g(s) > 1-7 



g(s)+M6 ~ I-7/2 
for all s £ [0,^(7)]. A direct computation shows that this is possible if we take 

75(0) 



S<8 { 







M(l -7)' 



Under this assumption r(x) provides a super-solution for (|3.39p and thus a positive solution of this 
problem can be constructed by the standard iteration procedure. 

In order to show that the sequence 4> m (x) is bounded we use C{ x )> the minimal positive solution 
of (13.391) and we need the following analog of the stratification Lemma 13.51 

Lemma 3.9 Fix M > 0. Let qA(x) be solution of 

-Aq A + Au ■ Vq A = M in ft, (3.42) 
qA = on <9ft. 

Then, given any 5 > there exist Aq > and Eq > so that for all solutions of {3.^2 ) we have 
< (Za(x) < 5 m G £ = ft \ ft £ for all A > Aq and e < Eq. 
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We do not present the proof of this lemma as it is essentially contained in that of Lemma 13.51 

We choose M > so that Ag(£(x)) < M for all x £ £l e , where C( x ) is the minimal positive 
solution of (I3.39p . We may also take A > sufficiently large, as in Lemma 13.91 We claim that then 
we will have, for all m > 1, (i) < 4> m (x) < 5 in G £ = 0, \ £l £ , and (ii) < 4> m (x) < C{x) for all 
x G Q £ . 

Let us prove this by induction. The function 4>i(x) satisfies 

-A0i + Au ■ V^i = Xg(0) in (3.43) 
01 = on 30. 

Our choice of M ensures that the right side in f)3.43|) is bounded above by M. Thus, if the amplitude 
A is sufficiently large we have < (f>i(x) < 5 in the tube G £ and in particular on d£l £ . Therefore, 
the difference si(x) = ((x) — <f>x(x) satisfies 

-Asi + Au ■ Vsi = X[g(qj) - g(0)] > in Q £ , (3.44) 
s\ > on dQ £ . 

It follows that si(x) > and < (f>x(x) < <^(x) in £ so that our claim holds for n = 1. Assume now 
that (i) and (ii) hold for n _i(x). It follows from the induction assumption that 

A 5 (0 n _i(x)) < \g(((x)) < M in O e , 

and A<7(0 n _i(x)) < A<7((5) < M in G £ . Therefore, for (f) n (x) we have 

-A0 n + Au ■ V(f> n = Xg(cj)n-i) < M, in O (3.45) 
4> n = on (90. 

Hence, by Lemma 13.91 we have < <j) m (x) < 5 in G £ . On the other hand, the difference s m (x) = 
((x) — (j) m (x) inside d£l £ obeys 

-As m + Au ■ Vs m = X[g(C) - g{<p m -i)] > in Q £ , (3.46) 
s m > on 90 e , 

and thus s m (x) > 0, so that < 4> m (x) < ((x) in Q £ . Therefore, the sequence 4> m is increasing and 
uniformly bounded from above. The limit 4> is a positive solution of f|3. 12j) and thus A < A* (^4). 
This finishes the proof of Lemma 13.81 □ 

3.6 General two-dimensional cellular flows 

We now look at the explosion problem 

-A(f) A + Au • Vc/) A = Xg{(/) A ) in O, (3.47) 

(p A = o on an, 

in a two-dimensional domain !JcK 2 with a cellular flow u which may contain more than one cell 
in 0, and complete the proof of Theorem 13.41 Let A^(^4) be the explosion threshold for the problem 
inside each cell: 

-A0 + Au ■ V0 = A#(0) in C,-, (3.48) 
(f) A = on Cj. 
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We already know that 



lim X*(A) = A*. 



from Proposition 13.61 and, of course, Xj(A) > X*(A) for all j. Hence, all we need to verify for the 
proof of Theorem [33] is that for any A < lim^^+oo Xj(A) solution of the problem (|3.47j) on the whole 
domain f2 exists. 

The proof is quite similar to the last part of the proof of Proposition 13.61 we construct the 
solution of (I3.47|) by the iteration procedure. Set 



Ao = min 

j 



lim X*(A) 



and take 7 > fixed. Consider any AG (0, (1 — 7)Ao), start the iteration process with 4>q = and 
define (j) m as the solution of 

-A<pi + An • V< = Xg{$t-i) m (3.49) 
6± = on dn. 



Lemma 3.10 The sequence <f>m( x ) is increasing in m, pointwise in x, for each A > Aq and has a 
uniformly bounded limit 4>a(x) £ L°°(f2) which is a solution of ftS^ty - 

The proof of Lemma 13.101 

The proof is quite analogous to that of Lemma 13.81 except that we use Lemma 13.51 where Lemma 13.91 
was used before. The first increment r]i(x) = (f>i(x) — 4>q(x) = (fri(x) satisfies 

-Ar/i + Au ■ V??i = As(0) > in ft, (3.50) 
771 = on (90. 

Hence, we have rji(x) > in and 4>i(x) > 4>o(x). Let us assume that r] n (x) = <fin( x ) ~ ( / ) n-i( x ) — 
in f2. As the nonlinearity g(s) is increasing in s, the function rj n+ i(x) is the solution of 

-Ari n+ i + Au ■ Vr] n+ i = X[g((f> n ) - g(<p n -i)} > in 0, (3.51) 
7]i = on d£l. 

Therefore, we have r] n+ i(x) > inside and thus the sequence (f>n( x ) ^ s increasing in n. We need 
to show that it is uniformly bounded from above. As in the proof of Lemma |3.8| we can find 5 > 
so that solution of the following explosion problem exists on each cell Cj for all A > Aq, 5 £ (0, 5q) 
and A < (1 — 7)Ao: 

— Aqj + Au • Vqj = Xg(qj) in Cj (3.52) 
qj = 5 on dCj, 

and, moreover, < qj(x) < -^(7) with the constant ^(7) which does not depend on the flow 
amplitude A. This is shown by exactly the same argument we used to show the existence of a 
positive solution for (|3,39|) . 

As qj(x) are uniformly bounded by ^(7), we may choose M > so that Xg(qj(x)) < M for all j 
and all We also take A > sufficiently large, as in Lemma [3. 51 (but with the right side of (|3.9p 

replaced by the constant M rather than 1). We claim that the following bounds will be preserved by 
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the iteration procedure: for all m > 1, (i) < 4> m ( x ) < 5 on the skeleton of separatrices T>o, and (ii) 
< 4> m ( x ) — Qj( x ) f° r an x ^ Cj- Again, we check this by induction. The function (j>i(x) satisfies 

-A#* + Au • V#* = Xg(0), in 0, (3.53) 

4>f = o on an. 

The right side in (|3.53j) is bounded above by M and hence, according to Lemma [3. 51 for the amplitude 
A sufficiently large we have < 4>i(x) < S on the skeleton T>q. Therefore, the first increment 
s\{x) = qj(x) — <f>\{x) satisfies 

-Asi + Au ■ Vsi = A[#(gj) - 5(0)] > in Cj, (3.54) 
si > on <9Cj, 

and thus si(x) > and < <j>\{x) < ^(2;) in Cj. Let us now assume that (i) and (ii) are true 
for 4> m _i(x) and show that they hold for 4> m (x) - the argument is exactly as for m = 1. By the 
induction assumption we have 

Xg(4> m -i(x)) < Xg(qj(x)) < M in the cell Cj, 

and thus 4> m (x) satisfies 

-A^ + An • V<pi = Xg(<pi^) < Xg{ qj (x)) < M, in Q, (3.55) 
^ = on da 

Hence, by Lemma [331 we have < (j) £ m {x) < 5 on "Do- Now, the difference s m (x) = 9j(ar) — (t> m { x ) 
inside Cj satisfies 

-As m + Ait • Vs m = X[g( qj ) - s(C-i)] ^ in C v ( 3 - 56 ) 
s m > on dCj, 

and thus s m (x) > 0, so that < (fi^ (x) < g'j(x) in Cj. Therefore, the sequence ^ is increasing and 
uniformly bounded from above. The limit <j> A (x) is a positive solution of (|3.47p . This completes the 
proof of Lemma 13.101 and thus that of Theorem 13.41 □ 
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